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Abstract. We resolve the paradox that although mag- 
netic collimation of an isotropic solar wind results in an 
enhancement of its proton flux along the polar directions, 
several observations indicate a wind proton flux peaked at 
the equator. To that goal, we solve the full set of the time- 
dependent MHD equations describing the axisymmetric 
outflow of plasma from the magnetized and rotating Sun, 
either in its present form of the solar wind, or, in its ear- 
lier form of a protosolar wind. Special attention is directed 
towards the collimation properties of the solar outflow at 
large heliocentric distances. For the present day solar wind 
it is found that the poloidal streamlines and fieldlines are 
only slightly focused toward the solar poles. However, even 
such a modest compression of the flow by the azimuthal 
magnetic field would lead to an increase of the mass flux 
at the polar axis by about 20% at 1 AU, relatively to 
its value at the equator, for an initially isotropic at the 
base wind, contrary to older and recent (Prognoz, Ulysses, 
SOHO) observations. For the anisotropic in heliolatitude 
wind with parameters at the base inferred from in situ ob- 
servations by ULYSSES/SWOOPS and SOHO/CDS the 
effect of collimation is almost totally compensated by the 
initial velocity and density anisotropy of the wind. This 
effect should be taken into account in the interpretation 
of the recent SOHO observations by the SWAN instru- 
ment. Similar simulations have been performed for a five- 
and ten-fold increase of the solar angular velocity corre- 
sponding presumably to the wind of an earlier phase of 
our Sun. For such conditions it is found that for initially 
radial streamlines, the azimuthal magnetic field created 
by the fast rotation focus them toward the rotation axis 
and forms a tightly collimated jet. 

Key words: Magnetohydrodynamics (MHD) - plas- 
mas - Sun: solar wind - stars: pre-main sequence - stars: 
winds, outflows - ISM: jets and outflows 
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1. Introduction 



Several stellar and extragalactic astrophysical systems 
have been observed to exhibit collimated outflows in the 
form of jets (young stellar objects, low and high mass X- 
ray binaries, black hole X-ray transients, symbiotic stars, 
planetary nebulae nuclei, supersoft X-ray sources, active 
galactic nuclei and quasars). In recent reviews of observa- 
tions from all such classes of astrophysical objects it has 
been argued that an interconnecting element may be a ro- 
tating accretion disk threaded by a magnetic field (Livio 
1999, Konigl & Pudritz, 1999). Such connection between 
the disk and the jet is particularly evident in HST obser- 
vations of several young stars in the nearby Orion nebula 
(Ray 1996). This rather convincing observational evidence 
of a close jet-disk relation is the basis for the presently 
prevailing view that an accretion disk is the necessary in- 
gredient for the production of collimated jets. 

On the other hand, theoretically it has been shown for 
quite some time by now that gas outflows from a rotat- 
ing magnetized object of any nature can be magnetically 
self-collimated to form a jet (Heyvaerts & Norman 1989, 
Chiueh et al 1991, Sauty & Tsinganos 1994, Bogovalov 
1995). This result seems to be a rather intrinsic property 
of magnetized winds with polytropic thermodynamics or 
not, where the self-compression of the plasma is provided 
by the toroidal magnetic field induced by the rotation 
of the central source. Henceforth emerges the generally 
accepted opinion that all observed jets are magnetically 
collimated (Livio 1999). However, no direct observational 
evidence exists today that most observed jets are indeed 
collimated solely by magnetic fields. Recently, it has been 
pointed out that the toroidal magnetic field is unstable 
and cannot collimate the jet effectively (Spruit et al. 1997, 
Lucek & Bell 1997) and it has been argued that magne- 
tized winds do not collimate without an external help, 
such as the channelling effects of a thick accretion disk 
and/or confinement from the ambient medium (Okamoto 
1999). It is thus crucially important to find direct observa- 
tional evidence that the magnetic field mainly collimates 
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the plasma in observed jets and by this way to test the 
theory of magnetic collimation. 

Nevertheless, plasma outflows do also emanate from 
isolated magnetized and rotating stars without an accre- 
tion disk, of which the solar wind (SW) is the classical 
and best studied example. The natural question which 
arises then is to what observable degree dynamical ef- 
fects are capable to collimate outflows from such single 
stars too. Theoretical studies on the angular momentum 
evolution of solar-type stars have concluded that at the 
end of the early accretion phase (PTTS) the star may be 
span up by more than 10 times the present solar rotation 
rate while its magnetic field is also strong (Bouvier et al. 
1997). And, in recent studies it has been shown that cold 
winds from such rapidly rotating and highly magnetized 
stars lead to considerable collimation of the outflow (Bo- 
govalov & Tsinganos 1999, henceforth Paper I). Similar is 
the result from studies of hot plasma outflows from effi- 
cient magnetic rotators (Sauty & Tsinganos 1994, Sauty 
et al. 1999). Hence, observation of the collimation effect 
in outflows from single stars could be the most reliable 
observational test of the theory of magnetic collimation. 

The question of the degree of collimation of the SW 
is an interesting possibility that has not been fully an- 
swered theoretically and observationally for quite some 
time now. Suess (1972) and Nerney & Suess (1975) were 
the first to model the axisymmetric interaction of mag- 
netic fields with rotation in stellar winds by a linearisation 
of the MHD equations in inverse Rossby numbers and to 
find a poleward deflection of the streamlines of the solar 
wind caused by the toroidal magnetic field. Later, Sakurai 
(1985) addressed the same problem by numerically solving 
the system of the polytropic MHD equations for the sta- 
tionary outflow. Washimi & Shibata (1993) modelled time 
dependent axisymmetric thermo-centrifugal winds with a 
dipole magnetic flux distribution on the stellar surface and 
a radial field in Washimi & Sakurai (1993). Polytropic 
MHD simulations of magnetized winds containing both a 
"wind" and a "dead" zone (Tsinganos & Low 1989, Mestel 
1999) have also been performed up to distances of 0.25 AU 
(Kcppcns & Goedbloed 1999). All these studies show some 
small deflection of the flow toward the axis of rotation. 

In the observational front, information on the degree 
of collimation of the SW can be inferred from anisotropics 
in the Lyman alpha emission. These solar UV photons 
are scattered by neutral H atoms of interstellar origin and 
where the SW mass flux is increased the neutral H atoms 
are destroyed and thus the Lyman alpha emission is re- 
duced. Early observations by the Mariner 10 (Kumar & 
Broadfoot 1979) and Prognoz (Bertaux et al. 1985) satel- 
lites have shown that there is less Lyman alpha emission 
near the equator in comparison to the ecliptic poles, than 
predicted by an isotropic SW (Bertaux et al. 1997). There- 
fore, these Lyman alpha observations imply that the SW 
mass efflux should be maximum at the equator and min- 
imum at the poles. The same trend is confirmed by in 



situ observations of Ulysses (Goldstein et al. 1996) and 
the SWAN instrument onboard of the SOHO spacecraft 
(Kyrola et al. 1998). However, the effect of SW collima- 
tion around the ecliptic poles would cause the opposite 
effect on Lyman alpha observations. In other words, al- 
though UV observations infer a SW mass efflux peaked 
at the equator, magnetic collimation would cause a SW 
mass efflux peaked at the poles, for an isotropic at the 
base wind. 

One of the main purposes of this paper is to resolve 
this paradox. We shall follow the idea of magnetic collima- 
tion of the SW and show which values of the parameters 
characterizing the heliolatitudinal dependance of the SW 
(Lima et al. 1997, Gallagher et al. 1999), such as density, 
bulk flow speed, mass efflux, etc are consistent with the 
observations by the Lyman alpha method. Furthermore, 
we shall follow the increase of the degree of collimation of 
a hot stellar wind by increasing the rotation rate of the 
star and show that a ten-fold increase of angular velocity, 
as is the case in the majority of the young rapid rotators, 
leads to a dramatic increase of the degree of stellar wind 
collimation. 

The paper is organised as follows. In Sect. 2 we justify 
the use of a split monopole model in our analysis for the 
collimation properties of the realistic solar wind at large 
distances from the Sun. In order to establish notation in 
Sect. 3 we give briefly the basic equations describing a sta- 
tionary and polytropic Parker wind. In Sect. 4, the initial 
configuration used together with the boundary conditions 
for the numerical simulation in the nearest zone are dis- 
cussed. In Sect. 5 the analytical method for extending the 
integration to unlimited large distances outside the near 
zone is briefly described while in Sect. 6 the parametriza- 
tion of the presented solutions is given. In Sect. 7 we dis- 
cuss the results for the isotropic SW in the near zone con- 
taining the critical surfaces and in the asymptotic regime 
of the collimated outflow, for a uniform rotation. In Sect. 
8 the case of a stellar wind from a star rotating faster than 
the Sun is taken up. A brief summary with a discussion 
of the main results is finally given in Sect. 9. 

2. The assumption of a split-monopole model for 
a stellar wind 

Fig. la is a sketch of the magnetic field structure in the 
corona of a star. Close to the stellar base the structure 
of the magnetosphere may be rather complicated. In this 
paper we are interested in the study of winds at distances 
much larger than the radius of the star where no closed 
field lines exist. It is therefore reasonable to consider the 
plasma flow starting at distances shown in Fig. la by the 
dashed line. For the solar wind the location of this surface 
can be put somewhere between the slow magnetosonic sur- 
face and the Alfvenic surface. We choose this location of 
the starting surface to avoid the solution of the problem of 
the wind acceleration in the very vicinity of the star which 
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Fig. 1. Schematic drawing of a white light image of the solar 
atmosphere indicating the shape of the solar magnetic field- 
lines. Beyond the dotted spherical distance located around the 
slow magnetosonic critical surface the magnetic fieldlines are 
radial to a first approximation. Electric sheet currents are lo- 
cated on the dotted lines. 

is defined not only by thermal pressure gradients but also 
by nonthermal processes of acceleration where the accel- 
eration mechanisms have not still studied sufficienty well 
and are beyond the scope of the present study (e.g., see 
Holzer & Leer 1997, Hansteen et al. 1997, Wang et al. 
1998). Above this base surface we can assume that the 
dynamics of the wind is mainly controlled by thermal and 
electromagnetic forces. In this approach the density and 
velocity of the plasma, together with the tangential com- 
ponents of the electric field and the normal component of 
the magnetic field are specified on this base surface while 
the tangential components of the magnetic field are free. 




(a) 



Fig. 2. Plasma flow from an axisymmetric rotator with an ini- 
tially split-monopole magnetic field, as in panel (a). Reversing 
the direction of the poloidal magnetic field in an arbitrary flux 
tube does not change the dynamics of the problem while we 
obtain the configuration shown in panel (b) which describes 
a nonstationary and nonaxisymmetric plasma flow from a ro- 
tator with a magnetic spot of opposite polarity on the base 
surface. The distribution of such spots can be arbitrary. 

Nevertheless, the solution of this problem is still too 
complicated. Open poloidal magnetic field lines go to in- 



finity and change their direction on the so called current 
sheets, some of which are indicated with dotted lines in 
Fig. la. These current sheets are present in any realistic 
wind from a stellar atmosphere, since the total magnetic 
flux of the open poloidal magnetic field lines is equal to 
zero while the mass loss rate is finite. The invariance prin- 
ciple summarized in the Appendix in a form appropriate 
to hot winds in a gravitational field allows us to simplify 
the structure of the magnetic field in the wind (Bogo- 
valov 1999). According to this principle, we can reverse 
the direction of some field lines so that the magnetic field 
is unipolar in each hemisphere, e.g., outward in the up- 
per and inwards in the lower hemisphere. In ideal MHD 
wherein we neglect all dissipative processes such as mag- 
netic reconnection, this operation does not affect the dy- 
namics of the plasma as long as the streamlines are not 
modified. This results in the configuration shown in Fig. 
lb where since we are not interested in the region up- 
stream of the base surface, this region is not shown. In 
this structure the current sheet is located only around the 
equatorial plane. To proceed, we further assume that the 
distribution of the normal component of the poloidal mag- 
netic field is uniform in the upper and lower hemispheres. 
In that case we get the model of the axisymmetric rotator 
with an initially split-monopole magnetic field. The field 
lines are magnetically focused toward the system's axis, 
as shown in Fig. 2a. 

We would like to stress that the model of the axisym- 
metric rotator describes not only axisymmetric outflows 
but also a wide class of nonstationary and axially non- 
symmetric outflows. This is due to the fact that according 
to the invariance principle (c.f. Appendix) the change of 
the direction of a magnetic field line in some flux tube 
docs not affect the dynamics of the plasma. For example, 
let's assume that we obtain a solution for the axisym- 
metric rotator, as shown schematically in Fig. 2a. Then, 
a reversal of the sign of some magnetic field lines in an 
arbitrary poloidal flux tube gives us a solution which is 
not axisymmetric and nonstationary, as shown in Fig. 2b. 
This is a solution for the plasma outflow from a rotator 
with uniform magnetic field at the base surface but with a 
magnetic spot of the opposite polarity on the upper hemi- 
sphere. Fig. 2b shows the cross-section of such a magnetic 
field by the poloidal plane. The stream lines arc the same 
as for the axisymmetric case. But the poloidal magnetic 
field changes sign in magnetic spots corresponding to the 
flux tube of the opposite polarity. The path of the field line 
in this flux tube in 3D is shown by a dashed line. These 
spots propagate in the poloidal plane with the velocity of 
the plasma and hence the pattern is nonstationary. 

It is clear that the number of such magnetic spots and 
their position at the base surface can be arbitrary. There- 
fore the study of the plasma outflow in the model of the ax- 
isymmetric rotator with an initially split-monopole mag- 
netic field allows us to study a much more wider classes 
of nonstationary and nonaxisymmetric flows. 
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3. The stationary polytropic Parker wind 

A Parker wind is taken as the initial state (t=0) for the so- 
lution of the time-dependent problem. In this initial state, 
the wind is assumed to flow along the radial magnetic field 
lines of an isotropic magnetic field (Parker 1963), although 
in general, the flow is not isotropic such that the wind has 
its own integrals of motion on every stream line ip = const. 
For simplicity, a polytropic relationship between the pres- 
sure P and the density p is assumed 



p = QWp 1 , 



where 7 is the polytropic index. Then, the Bernoulli equa- 
tion for energy conservation along a radial line in such an 
anisotropic wind from a nonrotating star has the form 
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2 + 7- 1 



,7-1 



GM 



(2) 



Denote by Voo(tp) the terminal velocity of the plasma 
on each fieldline. In order to get equations in dimensionless 
variables, we shall use for the radial distance R = R/R StCqi 
the density p = p/p Sieq , the entropy function Q = Q/Q eq 
and the velocity v = V/Vs,eq> in terms of the equatorial 
values of the sonic distance R s , C q, density p s ,cq, entropy 
function Q eq and sound speed V^ eq . 

The mass flux conservation in these dimensionless vari- 
ables takes the form 



pvR 2 = m(ip) , 

while the Bernoulli integral becomes, 



v 

T 



7-1 



QW>)p 
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GM 
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Since t | o(V') can be regarded as a function of R and p 
along a particular streamline ip, by taking the partial 
derivative of v 2 ^ with respect to p and R we obtain the 
usual Parker criticality relations which give the sound 
speed v s (tp) and the spherical distance Rs(ip) of the sonic 
critical surface along each streamline tp =const., 



v 2 s ^)= 1 Q^)p^ 1 \^). 
and 

i?m GM 1 



2V 2 cq R s , cq vlW) 



(5) 



(6) 



with the lower index s refering to the respective value of 
the variable at the sonic surface. Since, i? s ,eq = v s . cq = 1 
we have from the two criticality conditions, 



GM 



2^s 2 oq^s,eq 

such that 



= 1. 



v 2 (ip) 



(8) 



We are interested in obtaining the flow at large dis- 
tances from the central source. Therefore we shall take 
the distribution of the velocity and mass flux at infinity 
as the input parameters of the problem and introduce the 
parameter 



Taking into account the above two criticality conditions 
we have, 



(i) vtw^vZM 
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which gives 



and 
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5-37 
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(9) 



(10) 



(11) 



Note that Ucx>,eq > 1 only if 7 < 3/2. The enthalpy func- 
tion Q(ip) can be calculated in terms of £(ip) and m(ip) 
from Eq. (3) evaluated at the sonic surface and Eqs. (8) - 

(5), 



e(i>)_ = e^w 



^7-1 



(12) 



From Eqs. (3) - (8), the density at the critical surface is 
given in terms of £(-(/>) and m(?/>), 



PsW=e^)mW- 



(13) 



The Bernoulli equation in dimensionless variables has 
the form 



v 2 £ 5 " 37 W) 



1 



(l- 7 ) 



2 
R 



5-3 7 g 2 W 
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The above Bernoulli equation determines the plasma 
flow v(R;ip) along the prescribed radial magnetic field, 
once the polytropic index 7 and the distribution of the 
asymptotic velocity £(^>) are given. Then, Eq. (3) gives 
the density p = m/vR 2 once the mass flux m(ip) accross 
the poloidal streamlines is given. Note that in order to fi- 
nally calculate the physical variables V and p we need in 
addition, as input parameter of the problem, the equato- 
rial sound speed V s ^ eq while p SjCq can be calculated from 
the given m cq . 

Finally, consider the initial radial magnetic field 



(7) B = B, 



s,cq 



R 



(15) 



where B s cq is the magnetic field at the equatorial sonic 
transition. To define a dimensionless magnetic field, B, we 
need to normalize B to some characteristic value B c which 
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we choose to be given by the condition B 2 = 47rp Sjeq ^ eq . 
The dimensionless magnetic field B = B/B c then has the 
form 



density in the plasma is P 
the flow. 



Q(il>)p 7 , at any point along 



B = 
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(16) 



where VA,s,oq = -B SiCq /^/47rp s . oq is the AlfVenic velocity at 
the equatorial sonic point. Evidently, the strength of the 
initial dimensionless magnetic field is controlled by the 
magnitude of the ratio of the Alfven and sound speeds at 
the equatorial sonic distance, VA,s,eq/14. oq . 

4. The time-dependent stellar wind problem 

To obtain a stationary solution of the problem in the near- 
est zone of the star containing the sonic critical surface, 
it is needed to solve the complete system of the time- 
dependent MHD equations and look for an asymptotic 
stationary state. Then, the plasma flow in a gravitational 
field with the thermal pressure included is described by 
the following set of the familiar MHD equations, 



B p = 
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dt 
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(23) 



where we have used cylindrical coordinates (z, r, ip) and 
the magnetic field B has a poloidal magnetic flux denoted 
by ip(z,r). In the simulation we assumed a polytropic 
equation of state, so that the relationship of pressure and 



A correct solution of the problem requires a specifica- 
tion of the appropriate boundary conditions at the base 
surface. In this paper our main intention is to compare 
our results with observations at large distances from the 
source. In other words, we are interested in the asymptotic 
properties of the plasma flow in stellar winds. It is rea- 
sonable to start our numerical integration at a boundary 
surface R a placed just downstream of the slow mode crit- 
ical surface. In this way we avoid the problems connected 
with the correct description of the acceleration of the flow 
below the slow mode critical surface. Nevertheless, our 
boundary sourface will be placed below the Alfven and fast 
mode critical surfaces. The correct passage of the physi- 
cal solution from these two critical surfaces will yield the 
appropriate values of the toroidal component of the mag- 
netic field which is important in controlling the process of 
outflow collimation. 



Then, the appropriate physical conditions at the 
boundary surface of integration are in dimensional vari- 
ables: 

1. The density of the plasma at R — R a kept constant 
in time, although it may depend on the colatitude 9, 
p = p {6). 

2. The total plasma speed V o (0) in the corotating frame 
of reference at R = R , kept constant in time, although 
it may also depend on the colatitude 8, V 2 r ^ + V 2 Z s + 
(V(<p,o) ~ ^ r o) 2 = V D 2 (6') where V o (0) is the plasma speed 
on the surface R of the initial flow. The value of the initial 
velocity V {&) was taken such as to yield the observable 
values of the terminal velocity of the wind from a nonro- 
tating star, i.e., typical solar wind velocities at 1 AU. 

3. A constant and uniform in latitude distribution of the 
magnetic flux function at R = R Q , ip = ip . 

4. Finally, the continuity of the tangential component of 
the electric field across the stellar surface in the corotat- 
ing frame gives the last condition, (V( V;0 ) — 0.r o )B( p o ) — 
V(p,o)B( v , ) = 0. 



Recently it was realized by Ustyugova et al. (1999) that 
the boundary conditions at the outer box of simulation 
are important for obtaining the correct physical station- 
ary solution. The importance of a correct specification of 
the outer boundary conditions in MHD outflows has been 
emphasized previously in Bogovalov (1996, 1997) where 
we controlled the position of the outer boundary in the 
region where no signal can propagate from this boundary 
into the box of the simulation. For the solution of the full 
system of equations on the outer boundary we used only 
internal derivatives. 
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5. Method of numerical solution of the problem of 
stationary plasma flow to large distances from 
the star 

The asymptotic solution of the time-dependent problem 
in the nearest zone containing the Alfven and fast mode 
surfaces was used as the boundary condition for the far 
zone wherein we have a supersonic stationary flow. This 
boundary condition was then used in order to solve the 
system of the MHD equations describing the stationary 
outflow of supcrfast magnetosonic plasma. This system of 
equations consists of the set of the MHD integrals and of 
the transfield equation. 

As is well known, the stationary MHD equations ad- 
mit four integrals. They are the ratio of the poloidal mag- 
netic and mass fluxes, F(ip) = B p /4TrpV p ; the total angu- 
lar momentum per unit mass, rV v ~ FrB v — L(ip); the 
corotation frequency £l(ip) in the frozen-in MHD condi- 



tion V v B p - 



■V p B v = rB p Q(tp) and finally the total energy 



E(ip) in the equation for total energy conservation, 

The method of the solution of the transfield equation 
in the super fast magnetosonic region is described in de- 
tail in our previous Paper I. An orthogonal curvilinear 
coordinate system (ip, rj) is used, formed by the tangent 
to the poloidal magnetic field line fj = p and the first nor- 
mal towards the center of curvature of the poloidal lines, 
ip = V-0/|VV>|- A geometrical interval in these coordinates 
can be expressed as 



(dr) 2 = g 2 ,d^ 2 +g 2 d V 2 +r 2 d^ 2 



(24) 



the 



where g^,grj are the corresponding line elements, i.e., 
components of the metric tensor. 

According to Landau & Lifshitz (1975) the equation 
p-§${^jf)> where T tk is the energy-momentum 
tensor, will have the following form in these coordinates 
in the nonrclativistic limit, 



1 ip;k 
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(25) 



It is convenient to solve the transfield equation in the 
system of the curvilinear coordinates (tp, ry) introduced 
above. The unknown variables are z(r},tp) and r(r],tp). 
Therefore we need to know the quantities r v , z v , 

where r v = dr/drj, z n = dz/drj, = dr/dtp, z^ = dz/dtp 
and g n , g$. 

First, the metric coefficient g v is obtained from the 
above transfield equation (25), 



g v = exp 



(/ Gfa,V)#), 



(26) 



where 



G(r?,V) = 
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1 dr 
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The lower limit of the integration in (26) is chosen to 
be such that the coordinate rj is uniquely defined. In this 
way r\ coincides with the coordinate z where the surface 
of constant r\ crosses the axis of rotation. 

The metric coefficient g^ is given in terms of the mag- 
nitude of the poloidal magnetic field, 



9l> 



rB» 



(28) 



To obtain the expressions of r^p , z$, r n , z n we may use 
the orthogonality condition 



Tr/Tip ~t~ ZjjZip — , 



(29) 



and also the fact that they are related to the metric coef- 
ficients g v and g^, as follows, 
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2,2 
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(30) 



Finally, by combining the condition of orthogonality 



(29) and Eqs. (30) the remaining values of r^ , z Tj 
tained, 



ZipQrj 

9ip 



r\p9r, 
9ip 



are ob- 



(31) 



with g v above defined by the expression (26). For the nu- 
merical solution of the system of equations (31) the two 
step Lax-Wendroff method on the lattice with a dimension 
equal to 1000 is used. 

Eqs. (31) should be supplemented by appropriate 
boundary conditions on some initial surface of constant 
■q. The equations for r^, z^p defining this initial surface in 
cylindrical coordinates are, 



dr 

Ihp 



B z 
rB 2 > 



dz 
&4> 



(32) 



We need to specify on this surface the integrals 
E(tp), L(tp), Q(tp), F(ip) as the boundary conditions for the 
initial value problem. To specify the initial surface of con- 
stant j] and the above integrals, we use the results of the 
solution of the problem in the nearest zone when a station- 
ary solution is obtained for the time-dependent problem. 

6. Parametrization of the stationary solution 

It is convenient to consider the solution in dimcnsionless 
variables. In the present paper we express the velocities in 
units of the initial sound speed at the sound point on the 
equator V StCq , all the geometrical variables in units of the 
initial radius of the sound point on the equator i? s , cq and 
the magnetic field in units of B c = ^/47r/9 S! cqV^cq- 
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In this notation the solution depends on a few param- 
eters. Among them is the ratio of the radius of the star to 



the radius of the initial sound point i?* 
parameters 



R*/Rs 



s.cq 



and 



V a = 



V s . 



the 



(33) 



(34) 



together with the two dimcnsionlcss functions £("0) and 
fn{ip) which are equal to 1 in the case of uniform ejection 
of plasma at the base. 

Since we are not interested here in the dependence of 
the solution on the radius of the star, basically the flow is 
defined by the two parameters (3 and V a . 

Using the above definition of (3 and the relationships 
at the critical surface, Eqs. (7) - (8), this parameter can 
be expressed through observable parameters as 



QGM 5 

7-1 



37)3/2 



(35) 



where V^ is the terminal velocity of the plasma for the 
nonrotating star with mass M. The parameter 



V a = 



^Mi? s 2 eq V s , eq 



(36) 



where ip t = lim r ^oo B ■ r 2 is estimated on the equator of 
the nonrotating star and M is the total mass loss of the 
nonrotating star. 

In Paper I the degree of collimation of the outflow 
depended critically on a parameter a which was defined 
for cold plasmas as the ratio of the corotation speed at 
the Alfven transition to the initial velocity of the plasma. 
However for hot winds where the velocity depends on the 
radial distance we should introduce a more general defini- 
tion of this parameter a. 

In general, winds from astrophysical objects can be 
driven by a combination of several mechanisms of accel- 
eration such as the gradients of thermal pressure, Alfven 
wave and radiation pressures, magnetocentrifugal forces, 
etc. It is natural to characterize the efficiency of the 
magnetocentrifugal forces by the ratio V^/V^o where 
V m = (fi 2 V' t 2 /M) 1/3 is Michel's (1969) terminal velocity 
of a plasma accelerated only by magnetocentrifugal forces 
in the split-monopole model provided that the initial velo- 
city is zero while Voo. o is the terminal wind speed due to 
all other mechnisms of acceleration. In other words, Voo, o 
is the terminal velocity of the plasma for the nonrotating 
star. In this case the generalized a can be defined as 



' Vm 
Voofi 



3/2 



(37) 



In the special case of cold plasma outflow this defi- 
nition of the parameter a coincides with the parameter 
introduced in our previous Paper I. It can be easily shown 
that for polytropic winds this parameter can be expressed 
as 



a 



or, 



a 



7- 



1 



,3/4 



37' 



fan 



mv: 



3/2 



(38) 



(39) 



oo,0 



Magnetic rotators with a > 1 shall be called fast mag- 
netic rotators and those with a < 1 slow magnetic rota- 
tors. It is useful to note that for a slow magnetic rotator 
such as the Sun, R A « i/>t/(MV4o,o) 1/2 (MacGregor 1996), 
and hence a « Q.Ra/Voo, o- 

A quantitative classification of magnetic rotators on 
slow and fast (Belcher & MacGregor 1976) has also been 
introduced in Ferrcira (1997) by using the parameter 
CIRa/Va, where Va is the Alfven velocity in the Alfven 
transition located at the radius Ra- This parameter is also 
less than 1 for slow rotators, but for fast rotators it is of 
the order of 1, since for fast magnetic rotators Va ~ QRa 
(Michel 1969). Nevertheless, physically this classification 
of magnetic rotators to fast and slow coincides in both 
cases. Note also in passing that in terms of an energetic 
criterion for collimation deduced in Sauty & Tsinganos 
(1994) and Sauty et al. (1999) magnetic rotators are analo- 
gously classified as efficient (with cylindrical asymptotics) 
and inefficient (with radial asymptotics). 

7. Results for the isotropic and nonisotropic solar 
wind 

The most often used magnetized polytropic solar wind 
model is the classical one proposed by Weber & Davis 
(1967) where the poloidal magnetic field and stream lines 
are radial (see also MacGregor 1996). This model repro- 
duces the observed properties of the low speed streams 
at 1 AU within the observed fluctuations (Charbonneau 
1995). 

We will choose in our analysis the polytropic index 
to the value 7 = 1.1 such that the wind is heated. The 
spherical distance will be expressed in units of the distance 
of the slow magnctosonic point i? s i ow ~ 8.4 Rq and the 
velocity in units of the slow magnetosonic speed there, 
V s \ ow w 106 km/s. The fast magnetosonic transition occurs 
at i?f ast « 38 Rq where the fast magnetosonic speed is 
Vf as t ~ 230 km/sec. Note that for slow magnetic rotators 
like our Sun, the slow magnetosonic speed almost coincides 
with the sound speed and the Alfven critical point with 
the fast magnetosonic critical point. In particular, at the 
axis the Alfven and fast transitions coincide but at the 
equator the Alfven transition occurs earlier as in Figs. 3 
(see also Paper I and Kcppcns & Goedbloed 1999) 
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(a) 



G>) 



Fig. 3. Shape of poloidal magnetic field lines and streamlines in the near zone of an isotropic (panel a) and anisotropic (panel b) 
solar wind with 7 = 1.1, V a = 6.15 and f3 = 0.165. Spherical distance is given in units of radius of slow point at R s iow = 8.4 R & , 
with the base radius at R = 1.5 (dashed line). Thick lines indicate slow, Alfven and fast critical surfaces. 





Fig. 4. Shape of poloidal magnetic field lines in the far zone of an isotropic SW with 7 = 1.1 and f3 = 0.165. In (a) the poloidal 
field lines are plotted in a logarithmic scale, which magnifies their slight bending towards the axis, while in the linear scale of 
(b) it may be seen that collimation is negligible. 



In Fig. 3a we plot the shape of the poloidal magnetic 
field lines and streamlines in the near zone of a wind which 
is isotropic at the base with 7 = 1.1, V a = 6.15 and 
= 0.165. With these parameters a = 0.12, i.e., in our 
terminology the Sun is a slow rotator. Careful inspection 
of this figure shows that the flow is very slightly colli- 
mated to the axis of rotation. The solution in the far zone 
is presented in Fig. 4. 

In Fig. 4a the poloidal field lines of the SW are plotted 
in a logarithmic scale, which magnifies their slight bend- 



ing towards the axis. This logarithmic scale extends to 
the huge distance of 10 10 i? s i ow , i.e, about 4 x 10 8 AU « 60 
light years. This figure shows that the solar wind is indeed 
collimated toward the axis of rotation. But this collima- 
tion is indeed very weak. Fig. 4b shows the same field in a 
linear scale which shows that a jet is still not formed even 
at these huge distances. 

As is shown in Bogovalov (1995), the dependence of 
the magnitude of the poloidal magnetic field and den- 
sity on the cylindrical distance r becomes rather sim- 
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Fig. 5. Variation with dimensionless cylindrical distance X = r/R s \ ov of enclosed magnetic flux tp(X) (dot-dashed) and magnetic 
field B P (X) I B (solid) for f3 — 0.165 and 7 = 1.1 (panel a) and j3 = 0.044 and 7 = 1.2 (panel b). Dashed line indicates the 
analytically predicted values of B P (X)/B . 



pie if we assume for convenience that the MHD integrals 
E(ip), L(i/j), F(ip), and the terminal velocity of the 

jet Vj are constants and do not depend on ip and also that 
Vj S> Va(0), where V/t(0) is the Alfvenic velocity at the 
axis of rotation r = 0. In such a case, an approximate 
estimate of the dependence of the magnetic field on r is 
given in terms of the magnetic field and the density of the 
plasma on the axis of the jet, B p (0) and p(0), respectively, 



B p (r) _ p(r) _ 1 
B p (0) p(0) l + ir/Rjf 



(40) 



The radius of the core of the jet Rj is given in terms of 
the sound and Alfven speeds along the jet's axis V s (0) and 
Va(0), 



R i 




VaWJ n 



(41) 



Hence, the poloidal magnetic field and density remain 
practically constant up to axial distances of order Rj and 
then they decay fast like 1/r 2 outside the jet's core. 

The comparison of this theoretical prediction with the 
characteristics of the solar wind is shown in Fig. 5. There 
is a remarkable discreapancy between them which tells us 
that even at these unrealistically huge distance (the SW 
wind is presumably terminated earlier) the jet has still not 
formed. The reason is that the collimation of the wind oc- 
curs logarithmically with distance (Paper I) provided that 
the jet is not formed in the nearest zone, as it happens for 
fast rotators. For the solar parameters, there isn't enough 
radial distance to form the jet. 

In spite of the absence of a jet core in the solar wind, 
the lines of the plasma flow are certainly bent to the axis 
of rotation. And some observable effects can arise due to 



this bending. If the base density is isotropic, the SW mass 
efflux increases with the latitude 9 because of the magnetic 
focusing by about 20% from the equator to the pole (Fig. 
6a) at a distance from the Sun of about 5.8 AU. Approx- 
imately at this distance the interplanetary L a emission is 
formed, as observed by the SWAN instrument on board 
of SOHO. This theoretical anisotropy is in contradiction 
with the measurements of the anisotropy of the SW at 
large distance from the Sun by SWAN. This discreapency 
however, can be eliminated if we take into account some 
initial anisotropy in the SW. 

To study in more detail the effect of the focusing of the 
SW, we peformed calculations for a more realistic model 
of the SW including some initial anisotropy of the wind 
at its base. 

In Fig. 3b the shape of the streamlines and Alfven/fast 
critical surfaces is shown in the near zone of a wind which 
is latitudinally anisotropic in density and speed at the 
reference base at R = 1.5 R s , C q- For this simulation we 
used the latitudinal dependence of an exact MHD solution 
for the solar wind (Lima et al. 1997, Gallagher et al. 1999) 
with 



p(9)=p (l + S S m^8) , V r (e) = vJ ]^ Sm ll 6 (42) 

y 1 + 6 sin 9 

where Vo, and po correspond to the reference values of 
the flow speed and density and the distribution depends 
on the three parameters S, p and e. The values of e, 5 
and p have been chosen in Lima et al. (1997) as to best 
reproduce the observed values of the wind speed at various 
latitudes, as obtained recently by Ulysses (Goldstein et 
al. 1996). Their deduced values which we adopted in this 
study are 8 — 1.17, p = —0.38 and e = 8.6. A similar 
density enhancement about the ecliptic and an associated 
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(a) 



(b) 



Fig. 6. Distribution with latitude of characteristics of isotropic (panel a) and anisotropic (panel b) solar wind at 5.8 AU for 
(i = 0.165 and 7 = 1.1. Dashed lines indicate thermal pressure Pth(6)/Pth(0 = 0), solid lines mass flux J(9) = pvR 2 , dotted lines 
pressure of toroidal magnetic field Pt{0) / P t h(9 = 0) and dotted-dashed lines pressure of poloidal magnetic field P p (6) / P th (8 = 0). 



increase of the wind speed around the poles is also found 
in recent MHD simulations as well (Kcppens & Goedbloed 
1999). 

The distribution of the mass efflux at the distance of 
5.8 AU for the anisotropic at the base SW is shown in 
Fig. 6b. The effect of the magnetic focusing almost totally 
disappers at high latitudes. Near the equator the excess 
of the mass efflux remains remarkable in the region below 
30 degrees, although evidently the mass efflux decreases 
below 15 degrees. Therefore these results are in reasonable 
agreement with the distribution of the mass efflux of the 
SW which is deduced by SWAN at distances 5-7 AU. 

The drastical decrease of the effect of the focusing of 
the solar wind in the anisotropic case can be naturally 
explained by the larger velocity of the SW at high lat- 
itudes. According to observations by ULYSSES, the ve- 
locity at high latitudes is almost twice higher than the 
velocity at the equator (Feldman et al. 1996). It follows 
from the equation of motion Eq. (25) that the curvature 
radius R c of the streamlines in the poloidal plane at large 
distances can be estimated as (Paper I) 



1 8 r 2 Bl 



r 2 8tt rR c B p 



(43) 



The toroidal magnetic field at large distances can be esti- 



mated from the frozen in condition as B^ 
In this case we have 

1 

R 



8irpV 2 r 3 89 \ V P 1 p 



-(rQ/V p )B p . 



(44) 



The latitudinal distributon of pV p is almost constant 
with an ~ 15% increase only near the equator. Therefore 
the curvature radius depends on V p as 

1 1 (45) 



R. 



V 3 



provided that the mass flux density is fixed. Due to this 
strong dependence of the effect of collimation on the ve- 
locity of the plasma, the focusing practically disappears 
at high latitudes where the velocity is almost twice the 
velocity near the equator and the distribution of the mass 
efflux is in pretty good agreement with the observed one. 
The decrease of the mass efflux below 15° cannot be found 
by the present SWAN data analysis since it was assumed 
in this analysis that the mass efflux can only monotoni- 
cally increase with decreasing latitude. 

Up to now we have discussed the results of calcula- 
tions for the polytropic wind with 7 = 1.1. The flow with 
this polytropic index is almost isothermal and approxi- 
mates the SW near the Sun, up to several dozens of solar 
radii. But at larger distances the effective polytropic index 
should increase and at infinity it should go to the Parker 
value 3/2. At a first glance, it seems reasonable to ex- 
pect that the focusing of the plasma will be stronger for 
a higher polytropic index. 

To study this possibility we also performed calcula- 
tions for the isotropic solar wind with a polytropic index 
7 = 1.2, a reasonable value at the distance of several AU 
(Weber 1970). The change of the polytropic index results 
to a change of the other parameters f3 and V a so that the 
total magnetic flux, mass flux and terminal velocity of the 
plasma remain constant as in the case with 7 = 1.1. In 
particular, for 7 — 1.2 we have (3 — 0.044 and V a = 11.9. 
With these values of f3 and V a the parameter a remains 
constant. 

The distribution of the mass flux for 7 = 1.2 is shown 
in Fig. 7. 

It follows from this figure that the focusing of the flow 
to the axis of rotation becomes even smaller at distances 
below 5.8 AU than it does for 7 = 1.1. This is explained 
by the fact that the higher is the polytropic index, the less 
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Fig. 7. Distribution with latitude of the characteristics of the 
isotropic solar wind at 5.8 AU for /3 = 0.044 and 7 = 1.2. 
Dashed line indicates thermal pressure Pth(6) / Pth(d = 0), solid 
line mass flux J (9) — pvR 2 , dotted line pressure of toroidal 
magnetic field Pt(0)/Pth(9 = 0) and dotted-dashed line pres- 
sure of poloidal magnetic field P p (6) / Pth{9 = 0). 

is the variation of the plasma velocity with the distance. 
With a fixed terminal velocity, the wind with 7 = 1.2 has 
a higher velocity at the base than the wind with 7 = 1.1. 
But we already have discussed above the dependence of 
the focusing on the velocity of the flow. It follows from 
this discussion that the focusing is less for the wind with 
higher velocity. Therefore the focusing of the wind at sev- 
eral AU is maximal for 7 = 1.1. Since we have agreement 
of theory and observations of the solar wind anisotropy 
by SWAN for this polytropic index, it is expected that 
the wind with a higher polytropic index will not provide 
a higher level of collimation which could be inconsistent 
with the observations. 

The dependence of the focusing effect on the polytropic 
index becomes opposite at very large distances where the 
plasma velocity has already achived the terminal value. In 
this case it becomes important how fast the thermal pres- 
sure which tends to decollimate the plasma falls down with 
distance. The higher the polytropic index is, the faster the 
pressure falls down with the distance. Therefore at very 
large distances we should expect stronger collimation of 
the plasma for higher values of the polytropic index. This 
tendency is indeed found. Fig. 5 demostrates a comparison 
of the characteristics of the flow near the axis of rotation 
for polytropic indices 1.1 and 1.2. It is seen that the wind 
for 7 = 1.2 is stronger collimated although the formation 
of the jets is also not completed at these huge distances. 

8. Stellar winds from magnetic rotators faster than 
the Sun 

As we have seen in the previous section, the Sun is a rel- 
atively slow magnetic rotator with a < 1. A plausible 



senario for the Sun (and similar low mass stars) is that 
it has lost a large fraction of its angular momentum via 
magnetized a outflow while in its youth it was in a state 
of higher rotation, similar to the observed high rotation 
states of young stars (Bouvier et al. 1997) corresponding 
to larger values of a. It is interesting then to examine the 
change of the shape of the magnetic field of a stellar mag- 
netic rotator more efficient that the Sun, in the context 
of our simple modelling. For convenience and comparison 
with the present era Sun, we may keep constant the pa- 
rameters of the polytropic index and sound speed. In rapid 
rotators, the magnetic flux also increases roughly in pro- 
portionality to the rotation rate ft (Kawaler 1988, 1990). 
For simplicity let us neglect this increase and assume that 
the ratio of the Alfven and sound speed remains the same, 
V a = 6.15, as in the previous example, but the parameter 
(3 increases by 5 and 10 times, from j3 = 0.165 (a = 0.121) 
to /3 = 0.825 (a = 0.6) and then to = 1.65 (a = 1.21). 

We shall again consider in our numerical experiment 
that the star has a radial magnetic field and at t=0 rota- 
tion starts which via the Lorentz forces distorts this mag- 
netic field. After some time, a final equilibrium state is 
reached (Figs. 8) where the poloidal magnetic field and 
the plasma density are increased along the axis because 
of the focusing of the field lines towards the pole. A test 
that the steady state solution is reached is the constancy 
of the MHD integrals of motion, E(ip), L(ip), F(ip) and 
f2(V>)- Note that the wiggles appearing at the midlatitude 
in Figs. 8 are due to the step-like function representation 
of the surface of the star in the cylindrical coordinate sys- 
tem used in our calculations. This artifact is unavoidable 
in this system of coordinates and is also present in the 
calculations reported in Washimi & Shibata (1993). 

The shape of the poloidal field lines is shown in the 
near zone in Fig. 8. Collimation of the plasma to the axis 
of rotation already close to the source is evident. As in 
the cold plasma case of Paper I, the elongation of the 
subsonic region along the axis of rotation is present. The 
most surprising result is that this elongation occurs faster 
with an increase of the parameter a than it does for a 
cold plasma. It is easy to compare the flow in the nearest 
zone for a = 1.2 shown in Fig. 8b with the correspond- 
ing figure for similar a presented in Fig. 2 of Paper I. It 
appears that although the shape of the field lines is the 
same, the subsonic region for the cold plasma remains al- 
most spherical at this parameter in contract to the shape 
of the subsonic region for the hot plasma. The physics of 
this interest behaviour is as follows. 

The Alfven transition at a given point of the Z-axis oc- 
curs when V = (B p /A-KpV)B p . The ratio B p /pV remains 
constant along a field line. Therefore the right hand side 
in this expression increases with collimation as B p . In the 
cold plasma limit V is constant and the displacement of 
the Alfven point down the flow in the cold plasma occurs 
only due to an increase of B p . In the hot plasma case, the 
collimation also modifies the velocity V. 
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7=1.1 
/?=0.825 




(a) 



(b) 



Fig. 8. Shape of poloidal magnetic field lines/streamlines in near zone of isotropic wind from a star rotating 5 times (/3 = 0.825, 
panel a) and 10 times faster than the Sun, (/3 = 1.65, panel b) for 7 = 1.1 and V a = 6.15. Spherical distance is in units of radius 
of slow point at R s iow = 8.4 Rq, with the base radius at 1.57i s i ow (dashed line). The initial nonrotating monopole magnetic field 
has a spherical Alfven surface at 4.5R B i ov . Thick lines indicate slow, Alfven and fast critical surfaces. 
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Fig. 9. Flow speed at the axis (solid) and the equator (dashed 
line) for wind from a star rotating 10 times faster than the 
Sun, with 7 = 1.1 and f3 — 1.65. 

This modification is shown in Fig. 9. The solid line 
shows the velocity at the axis and the dashed line shows 
the velocity at the equator. Collimation results in an de- 
crease of the thermal pressure gradient. Therefore the 
thermal acceleration of the plasma becomes less efficient. 
The velocity of the plasma even slightly decreases due to 
the gravitation force. Therefore the Alfven surface in the 
hot case elongates along the axis of rotation with an in- 
crease of a, faster than it does in the cold plasma case. 
It is clear that at some value of a specific for every flow, 
the subsonic region near the axis will be elongated to in- 
finity in the Z-dircction. In the hot plasma this happens 
at smaller a than it does in the cold plasma case. This 



subfast region should be certainly unstable. Therefore, 
the plasma flow at these conditions cannot be stationary. 
Those jets with a ^> 1 eventually should become turbu- 
lent with properties differing from the properties of the 
stationary jets found for a < 1. The physics of these jets 
should be examined in a separate study. 

In Fig. 10 the shape of the poloidal ficldlines is shown 
for the case of the rapid rotator. Their bending towards 
the rotation axis is evident not only in the logarithmic plot 
of Fig. (10a) but also in the linear plot of Fig. 10b where 
a tightly collimated jet is formed already at a relatively 
small distance from the source. 

In Figs. 11 the poloidal component of the magnetic 
field is plotted together with the magnetic flux enclosed 
by a cylindrical distance X. The poloidal magnetic field 
B p (X)/B (solid line) is given in units of its reference 
value B 0l corresponding to the magnetic field at the sym- 
metry axis r = and some reference height Z a = 5 • 10 s for 
(3 = 0.825 (Fig. 11a) and Z a = 1.66 x 10 6 for (3 = 1.65 (Fig. 
lib). The asymptotic regime of the jet is achived at these 
distances. The dashed curve gives the analytically pre- 
dicted solution for the poloidal magnetic field B P (X)/B , 
Eq. (42). Note that the agreement between the calculated 
and analytically predicted values of the radius of the jet is 
pretty good. It follows from this comparison that for these 
parameters we indeed achieve the distance where the jets 
are really formed. 

9. Discussion of the results 

One of the basic results from the numerical simulations 
reported in this paper as well as in Paper I is that plasma 




Fig. 10. Shape of poloidal magnetic field lines for a wind from a star rotating 10 times faster than the Sun, with 7 = 1.1 and 
P — 1.65. In (a) the poloidal field lines are plotted in a logarithmic scale, which artificially magnifies their bending towards the 
axis, while in the linear scale of (b) a higher degree of collimation in comparison to the corresponding case of Fig. 4b may be 
seen. 




Fig. 11. Variation with dimensionless cylindrical distance X — r/R s i ov , of enclosed magnetic flux ip(X) (dot-dashed) and 
magnetic field B P (X) / B (solid) from a moderate stellar magnetic rotator with f3 = 0.825 (panel a) and a faster magnetic 
rotator with /3 = 1.65 (panel b). Dashed line indicates the analytically predicted values of B P (X) / B a . 



ejected by rotating magnetized objects has the property 
of self-collimation. No other special conditions are needed 
to get a collimated outflow. Therefore, jets should be a 
rather common phenomenon in astrophysics. This is the 
most important prediction of the theory of magnetic col- 
limation which is in pretty good agreement with observa- 
tions, since numerous jets are observed to be associated 
with astrophysical objects of different nature. 



However, a direct application of the theory of magnetic 
collimation to an isotropic solar wind predicts an increase 
of the mass efflux near the poles while observations of the 
distribution of the solar wind mass efflux with heliolati- 
tude at distances 5-7 AU by a range of instruments have 
given the opposite result. In this work we eliminated this 
apparent contradiction by taking into account in the calcu- 
lations more realistic parameters for the anisotropic solar 
wind, as those inferred recently by ULYSSES and SOHO. 
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Then, a recalculation of the heliolatitudinal distribution 
of the SW mass efflux at large distances with an initially 
anisotropic distribution of the density and flow speed at 
the base show that the initial excess of the mass flux near 
the equator is preserved up to the outer boundary of the 
SW. Magnetic collimation is remarkable only in the region 
of the lower speed wind near the equator. It results in a 
relative decrease of the mass efflux in comparison with the 
flux at the low latitudes less than 15°. Can this effect be 
found in the SWAN data? To answer this question a de- 
tailed comparison of the calculations with these data is 
necessary. 

It is widely believed that due to the observed close 
disk-jet connection collimated outflows are possible only 
from systems containing accretion disks. It is clearly 
demonstrated in this study that fast rotating stars without 
disks can also produce jet-like outflows. This prediction is 
crucially important for the theory of magnetic collima- 
tion, since it gives a most reliable observational test of the 
theory. Up to now there has been no direct observational 
evidence that observed jets are really collimated by the 
magnetic field. And, observation of jet-like outflows from 
objects which do not contain accretion disks would be this 
needed evidence. 

Young fast rotating stars of solar mass are among the 
objects which produce jets. The effect of magnetic self- 
collimation of the winds is parameterized by a, which can 
be presented in the form 



0.12 



(^/V> )(n/n Q ) 
yV'//M )(W^)) 3/2 ' 



(46) 



for Vq « 400 km/s. Some young stars with solar mass 
have angular velocities up to 100 times the angular velo- 
city of the Sun (Bouvier et al. 1997). For example, in AB 
Doradus ft* = 54 £1 Q (Jardine et al. 1999). A phenomeno- 
logical dynamo mechanism (MacGregor 1996) predicts a 
magnetic flux which scales as ip*/ipQ oc Q*/Qq. However, 
the dependence of the mass flux on is not known. The 
theoretical analysis in the Weber & Davis (1967) approxi- 
mation shows that the mass flux is practically independent 
of the angular velocity. In such a case AB Doradus would 
have a w 350. On the other hand, if the mass loss rate is 
proportional to the magnetic pressure, M*/M Q oc B^/Bq, 
a value of a which is 54 times smaller results, i.e., a « 6.5. 
This star is certainly a rapid rotator and should produce a 
collimated outflow. However, can such jets be observed ? 
Unfortunately their mass flux may be too small and thus 
it may be practically impossible to observe these jets di- 
rectly. 

Outflows from stars with much higher mass loss rates 
could be modelled with this study. For example, B/Be 
stars have massive radiation driven winds. Usually these 
stars also rotate rapidly, at least Be stars. Their mass 
loss rate lies in the range M* = (10~ 6 — 10~ 8 )M Q /year, 
while their angular velocity is about ft* ~ 20f2 Q . The 



average magnetic field on the surface of B/Be stars can 
vary from 200 G to 1600 G (MacGregor 1996). Assuming 
that i?* = lOi?©, the magnetic flux from such a star is 
V>* = (7xlO 3 -5.6xlO 4 )V'0. With wind speeds of the order 
of Voo = 1000 km/s, the parameter a lies in the range 
a = 0.05 — 5. This means that some of these B/Be stars 
(but not all) could be fast magnetic rotators producing 
therefore jet-like outflows. Due to their huge mass loss 
comparable to the mass loss of classical T Tauri stars, 
these jets could be more easily observable. It is interesting 
in this connection to note observations by Marti et al. 
(1993) of jets from a B star in the HH 80/81 complex. 

The most strikingly unwanted result obtained in this 
study is that the part of the total magnetic flux (and mass 
flux correspondingly) going in to the jet is only of the order 
of 1%, as it may be seen in Fig. 11. Almost all magnetic 
flux goes in to the radially expanding wind with mass loss 
M w . In this case the results of our calculations can not 
be directly applied to jets from YSO because in most of 
them it appears that the mass flux in the jets Mj is about 
1% of the accretion rate M a (Hartigan et al. 1995). If we 
assume that only about 1% of the outflowing wind goes 
in to the jet, like in our results, then we will have the 
uncomfortably high ratio M w /M a ~ 1, which apparently 
should not be so high for outflows from accretion disks 
(Pellcticr & Pudritz, 1992). This means that some of our 
assumptions are not valid in the central source of jets from 
YSO. Modifications of the input parameters of the model 
which provide collimation of an arbitrary high fraction of 
the wind into the jet will be discussed in our next paper. 

Finally we would like to emphasize once more the sit- 
uation concerning wind flows from sources with a 3> 1. It 
seems that the majority of interesting sources of outflows, 
such as rapidly rotating stars and systems with accretion 
disks are just those which satisfy this condition. Jets from 
such sources should be nonstationary and turbulent with 
properties which may strongly differ from the properties 
of the laminar jets conidered in this paper. The physics of 
such jets still remains to be studied. 
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A. Appendix 

We show that the dynamics of plasma in ideal MHD flows 
is invariant in relation to a reversal of the direction of the 
magnetic field lines in an arbitrary flux tube. Firstly this 
property of ideal MHD flows was used for the solution 
of the problem of plasma outlow from oblique rotators in 
pulsar conditions (Bogovalov 1999). Here we simply show 
that it is also valid for flows in a gravitational field with 
thermal pressure. 
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The plasma flow in the nonrelativistic limit is de- 
scribed by the familiar set of the ideal MHD equations 

p ilt + p( - v ' V ^ V = ~ VP ~ pV$ + x B ^ x B '^ A1 ) 



<9B 



= V x (V x B), 



VB = 0, 



dp 
dt 



+ V • (pV) = 0, 



(A2) 
(A3) 



(A4) 



where $ is the gravitational potential and P is the pres- 
sure. 

Let us assume that we have some solution which is de- 
scribed by the functions B(r, t), p(r, t), V(r, t) and P(r, t). 
We show that the change of the direction of the magnetic 
field in an arbitrary magnetic flux tube does not change 
the dynamics of the plasma. 

Let us introduce a scalar function rj(r, t) with the prop- 
erty that rj = 1 everywhere except inside the choosen flux 
tube where r] = —1. This function satisfies the following 2 
conditions 



B • V?7 = , 



and 
drj 



+ V • Vr? = . 



(A5) 



(A6) 



The second equation is the consequence of the frozen-in 
condition such that the value of r\ is advected together 
with the plasma. 

Then the solution rjB{r,t), p(r,t), V(r,i) and P(r,t) 
also satisfies the system of Eqs. (A1-A4). Indeed, the 
Lorentz force in the right hand side of (Al) is, 

[V x (77B)] x (jjB) = 7?[V7? xB + r ? VxB]xB = 
V(r/ 2 /2) x B + ry 2 (V x B) x B = (V x B) x B , (A7) 

since rj 2 = 1. This means that the forces affecting the 
plasma do not change with this transformation. Eqs. (A2, 
A3) are satisfied due to conditions (A5, A6). 
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